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Measurements of HERA that explore the parton dynamics at low Bjorken x are 
presented together with precise determinations of the strong coupling constant a s . 
Calculations at next to leading order using the DGLAP evolution fail to describe 
the data at low x and forward jet pseudorapidities. The a s {Mz ) measurements at 
HERA are in agreement with the world average and have very competitive errors. 


1. Introduction 

At HERA, forward jet production and jet-jet correlations are expected to 
be sensitive to the parton dynamics at low Bjorken x (aiBj)- In the first 
section, these type of measurements and comparison to DGLAP, BFKL 
and other models are presented. For the measurements presented in the 
second section, jets were selected in kinematic regions where the proton 
parton distribution functions (PDFs) are well constrained and the DGLAP 
equations are valid. These measurements allow precise tests of perturbative 
QCD (pQCD) and the determination of the strong coupling constant, a s . 


2. QCD dynamics at low x 

Inclusive forward jet production was measured by the HI Collaboration 1 for 
events in the kinematic region 5 < Q 2 < 85 GeV 2 and 10 -4 < XBj < 4-10 -3 , 
where Q 2 is the virtuality of the exchanged photon. Jets were found in the 
laboratory frame with >3.5 GeV, 7° < 0j e t < 20° (corresponding to 
the 1.74 < < 2.8 a ), Xj e t > 0.035, and 0.5 < E^,- et /Q 2 < 5. E^ is the 

transverse energy of the jet and Xj e t is the fractional energy of the proton 
taken by the jet. 


a T) — — log[tan(0/2)] is the pseudorapidity, where 0 is the polar angle. 
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The measured cross section is shown in Fig. 1 as a function of Sej and 
compared to the prediction of NLO calculations from DISENT (left) and 
various QCD models (middle). The DISENT calculations were performed 
using the CTEQ6M parametrization of the proton PDFs. The renormal¬ 
ization scale was chosen to be = E F t . The DGLAP model with direct 
photon interactions alone (RG-DIR, RAPGAP) and the NLO calculation 
fall below the data, especially at low XBj- The description of the data by 
RAPGAP is significantly improved if contributions from resolved photon 
interactions are included (RG-DIR+RES). The Color Dipole Model (CDM) 
shows a similar behaviour to RG-DIR+RES. In addition, the CCFM based 
CASCADE MC predicts a different shape of the distribution that results 
in a poor description of the data. 



Figure 1. Cross section for inclusive forward jet production as a function of x com¬ 
pared to NLO calculation (left) and QCD MC models (middle). Measured differential 
cross section (right) for inclusive jet production in the forward jet rapidities as a function 
of X B j • 


Jet production in NC DIS has been measured by the ZEUS Collabora¬ 
tion 2 for Q 2 > 25 GeV 2 , y > 0.004, E' e > 10 GeV (where E' e is the energy 
of the scattered positron), and cos( 7 h) < 0 b . Jets were selected in the lab¬ 
oratory frame with E^ > 6 GeV, 2 < rf et < 3, and 0.5 < (E^) 2 /Q 2 < 2. 

The measurement is presented in Fig. 1 (right) as a function of iBj and 
compared to NLO calculations (DISENT, y. R = /j F = Q, CTEQ6) and 
to the CDM (ARIADNE) and MEPS (Matrix Elements + Parton Showers, 
LEPTO) models. The CDM prediction gives a reasonable description of the 


b 7 h is the hadronic angle. It corresponds, in the Quark Parton Model, to the angle of 
the scattered quark. 
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data. The NLO and MEPS predictions fail to describe the data in the low 
a:Bj region. The uncertainty induced by the variation of the renormalization 
scale is large, indicating that missing higher order or ln(l/:r) terms in the 
calculation could be important in this region. 

The HI Collaboration has measured inclusive dijet production in DIS 3 
in the kinematic range 5 < Q 2 < 100 GeV 2 , 10 -4 < x < 10 -2 , and 
0.1 < y < 0.7. Dijets were reconstructed in the hadronic center-of-mass 
system (HCM) and selected with the requirements: —1 < r?j e t,iab < 2.5 and 
(E^ n ’ 2 )* >7,5 GeV. The azimuthal asymmetry is defined by: S(a ) = 
ri8o° Vd ‘ ]et ( A '^ > X ^ ■ >dA< ^ > — where A<^>* is the azimuthal separation in the 

Jo A r dijet(A0*,a:,<2 2 )dA0* 

HCM frame between the two hardest transverse energy jets. 

Figure 2 presents the S distribution for a = 120° as a function of x for 
different Q 2 compared to predictions of DGLAP NLO calculations (left) and 
different models (right). The measured values of S for the chosen a are of 
the order of 5% and increase with decreasing x. The rise is most prominent 
in the lowest Q 2 bin. The DISENT (0(a 2 )) calculations predict no rise of 
S with x and fall below the measurements. The NLOJET (0(a 3 )) give a 
good description of the data at large Q 2 and large x but fail to describe 
the strong rise towards low x, particularly in the lowest Q 2 range. A good 
description of the S distribution at low x and low Q 2 is predicted from the 
CDM model. The CCFM based CASCADE predictions depend strongly 
on the choice of the unintegrated gluon density. 



Figure 2. Azimuthal asymmetry S for A<f>* < 120° as a function of ccbj and Q 2 . The 
data are compared to NLO predictions (left) and to the CCFM and CDM models (right). 
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3. Precise tests of QCD and the measurement of ot s 



Figure 3. Inclusive dijet and trijet cross sections (left) and their ratio as a function of 
Q 2 (middle) together with the extracted values of Q! S in each Q 2 bin. Mean integrated 
jet shape (ifj(r = 0.5)) versus and rf et (right) in comparison with the expectation 
for three different a s values. 


th. uncert. 
exp. uncert. 


Differential dijet and trijet cross sections have been measured by the 
ZEUS Collaboration 4 in the kinematic range 10 < Q 2 < 5000 GeV 2 and 
0.04 < y < 0.6. Jets were found in the Breit frame and dijet (trijet) events 
were selected with: —1 < r/j e t,iab < 2.5, Et, jet,Breit > 5 GeV, and invariant 
mass M 2 / 3 jets > 25 GeV. 

Figure 3 shows the measured di¬ 
jet and trijet cross sections as a func¬ 
tion of Q 2 (left) and their measured 
ratio i? 3 / 2 (middle). The NLO pre¬ 
dictions of NLO JET {fin = fip = 

{E 2 t + Q 2 )/4, CTEQ6) corrected for 
hadronization effects are compared to 
the data. NLOJET provides a good 
description of both the shape and the 
magnitude of the measured cross sec¬ 
tion. The correlated systematic and 
the renormalization scale uncertain¬ 
ties largely cancel in the ratio of the 
cross sections. This cancellation al¬ 
lows the extraction of a s {Mz) with a 

good precision down to Q 2 of 10 GeV 2 , using a method similar to that 
of a previous ZEUS publication 5 . The value of a s was measured to be 
oc s {Mz) = 0.1179 ± 0.0013 (stat.) 1“ (exp.) t°;“° 4 (th.). 


Jet shapes in NC DIS 

ZEUS (Nucl Phys B 700 (2004) 3) 

Multi-jets in NC DIS 
ZEUS (DESY 05-019 - hep-ex/0502007) 
Inclusive jet cross sections in vp 
ZEUS (Phys Lett B 560 (2003)7) 

Subjet multiplicity in CC DIS 
ZEUS (Eur Phys jour C 31 (2003) 149) 
Subjet multiplicity in NC DIS 
ZEUS (Phys Lett B 558 (2003) 41) 

NLO QCD fit 

HI (Eur Phys J C 21 (2001) 33) 

NLO QCD fit 

ZEUS prel. (contributed paper to ICHEP04) 
NLO QCD fit 

ZEUS (Phys Rev D 67 (2003) 012007) 
Inclusive jet cross sections in NC DIS 
HI (Eur Phys J C 19 (2001) 289) 

Inclusive jet cross sections in NC DIS 
ZEUS (Phys Lett B 547 (2002) 164) 

Dijet cross sections in NC DIS 
ZEUS (Phys Lett B 507 (2001) 70) 

World average 

(S. Bethke. hep-ex/0407021) 


Figure 4. A summary of the from 
measurements at HERA. 
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Jet substructure for jets produced in DIS has been measured by the 
ZEUS Collaboration 6 for Q 2 > 125 GeV 2 . Jets were reconstructed in the 
laboratory frame and were selected with UjA >17 GeV and — 1 < rf et < 
2.5. The mean integrated jet shape is defined as the averaged fraction of the 
jet transverse energy inside the cone of radius r: {ip(r)) = JY ets , 

where E ^ is the total transverse energy of the jet, Et(t) the part of it 
inside the cone of radius r and Nj eis is the total number of jets in the 
sample. 

Figure 3 (right) shows the measured {i/j(r = 0.5)) as a function of 
ryj et and E 3 ^. There is no significant variation of (ip(r = 0.5)) with if et 
in DIS, whereas (ip(r = 0.5)) increases as E^ increases. The sensitiv¬ 
ity of the measurements to the value of a s (Mz) is illustrated in Fig. 3 
(lower part of plot) by comparing the measured = 0.5)) to NLO 
QCD calculations using three different values of a s (Mz). The NLO QCD 
calculations provide a good description of the measured (ijj(r = 0.5)) 
and thus this observable was used to determine a s (Mz). The value of 
a s (Mz) as determined by fitting the NLO QCD calculations to the mea¬ 
sured mean integrated jet shape (ip(r = 0.5)) for E 3 ^ > 21 GeV is 
a s {M z ) =0.1176 ±0.0009 (stat.) (exp.) (th-)- 

4. Conclusions 

The HERA measurements show that DGLAP NLO calculations at low x 
and forward jet pseudorapidities fail to describe data but the large theo¬ 
retical uncertainties prevent a decisive conclusion on parton dynamics at 
low x. A summary of the a s measurements at HERA is shown in Fig. 4. 
The a s measurements at HERA are in agreement with the world average 
and have very competitive errors. For more accurate measurements of a s 
improved theoretical calculations would be needed. 


References 

1. HI Collab., Contributed Paper N 5-0172 to 32nd ICHEP, Beijing, 2004 

2. ZEUS Collab., DESY-05-017 (2005) 

3. HI Collab., A. Aktas et al., Eur. Phys. J. C37 141 (2004) 

4. ZEUS Collab., DESY-05-019 (2005) 

5. ZEUS Collab., Phys. Lett. B 507, 70 (2001) 

6. ZEUS Collab., Nucl. Phys. B 700, 3 (2004) 



